Predicting the future of our oceans-Evaluating genomic forecasting approaches in marine species.

Climate change is restructuring biodiversity on multiple scales and there is a pressing need to understand the downstream ecological and genomic consequences of this change. Recent advancements in the field of eco-evolutionary genomics have sought to include evolutionary processes in forecasting species' responses to climate change (e.g., genomic offset), but to date, much of this work has focused on terrestrial species. Coastal and offshore species, and the fisheries they support, may be even more vulnerable to climate change than their terrestrial counterparts, warranting a critical appraisal of these approaches in marine systems. First, we synthesize knowledge about the genomic basis of adaptation in marine species, and then we discuss the few examples where genomic forecasting has been applied in marine systems. Next, we identify the key challenges in validating genomic offset estimates in marine species, and we advocate for the inclusion of historical sampling data and hindcasting in the validation phase. Lastly, we describe a workflow to guide marine managers in incorporating these predictions into the decision-making process.

Differential expression and alternative splicing analyses of multiple tissues reveal albinism-associated genes in the Wels catfish (Silurus glanis).

Albinism is a widespread departure from a typical body colouration due to altered melanin production. The Wels catfish (Silurus glanis) is among the largest freshwater fish species in the world, and albino individuals occur both in the wild and in aquaculture. Here, we performed transcriptome-wide analysis of albino and normally pigmented S. glanis using four tissues (skin, dorsal fin, whole eye and liver) to identify genes associated with albinism by exploring patterns of differential expression (DE) and differential alternative splicing (DAS). Multi-tissue analyses revealed a large number of genes in skin (n = 1355) and fin (n = 614) tissue associated with the albino phenotype in S. glanis, while the number of DE genes in eye and liver tissues was lower (n = 188, n = 189, respectively). Several DE genes across multiple tissues were detected as the most promising candidates (e.g., hsp4, hsp90b1, raph1, uqcrfs1, adcy-family and wnt-family) potentially causally linked to the albino phenotype in Wels catfish. Moreover, our findings supported earlier observations of physiological differences between albino and normally pigmented individuals, particularly in energy metabolism and immune response. In contrast, there were only a few pigmentation-related genes observed among DAS genes (4 in skin, 2 in fin), the overlap between DAS and DE genes was low (n = 25) and did not include known pigmentation-related genes. This suggests that DAS and DE in Wels catfish are, to a large extent, independent processes, and the observed alternative splicing cases are probably not causally linked with albinism in S. glanis. This work provides the first transcriptome-wide multi-tissue insights into the albinism of Wels catfish and serves as a valuable resource for further understanding the genetic mechanisms of pigmentation in fish.

Distribution and prevalence of the myxozoan parasite Tetracapsuloides bryosalmonae in northernmost Europe: analysis of three salmonid species.

Global climate change is altering the abundance and spread of many aquatic parasites and pathogens. Proliferative kidney disease (PKD) of salmonids caused by the myxozoan Tetracapsuloides bryosalmonae is one such emerging disorder, and its impact is expected to increase with rising water temperature. Yet, the distribution and prevalence of T. bryosalmonae in Northern Europe remain poorly characterized. Here, we studied 43 locations in 27 rivers in northernmost Norway and Finland to describe T. bryosalmonae infection frequency and patterns in 1389 juvenile salmonids. T. bryosalmonae was discovered in 12 out of 27 rivers (44%) and prevalence ranged from 4.2 to 55.5% in Atlantic salmon and from 5.8 to 75% in brown trout among infected rivers. In sympatric populations, brown trout was more frequently infected with T. bryosalmonae than was salmon. Age-specific parasite prevalence patterns revealed that in contrast to lower latitudes, the infection of juvenile fish predominantly occurs during the second summer or later. Temperature monitoring over 2 yr indicated that the mean water temperature in June was 2.1 to 3.2°C higher in rivers containing T. bryosalmonae compared to parasite-free rivers, confirming the important role of temperature in parasite occurrence. Temporal comparison in T. bryosalmonae prevalence over a 10 yr period in 11 rivers did not reveal any signs of contemporary parasite spread to previously uninfected rivers. However, the wide distribution of T. bryosalmonae in rivers flowing to the Barents Sea indicates that climate change and heat waves may cause new disease outbreaks in northern regions.

Association mapping reveals candidate loci for resistance and anaemic response to an emerging temperature-driven parasitic disease in a wild salmonid fish.

Even though parasitic infections are often costly or deadly for the host, we know very little which genes influence parasite susceptibility and disease severity. Proliferative kidney disease is an emerging and, at elevated water temperatures, potentially deadly disease of salmonid fishes that is caused by the myxozoan parasite Tetracapsuloides bryosalmonae. By screening >7.6 K SNPs in 255 wild brown trout (Salmo trutta) and combining association mapping and Random Forest approaches, we identified several candidate genes for both the parasite resistance (inverse of relative parasite load; RPL) and the severe anaemic response to the parasite. The strongest RPL-associated SNP mapped to a noncoding region of the congeneric Atlantic salmon (S. salar) chromosome 10, whereas the second strongest RPL-associated SNP mapped to an intronic region of PRICKLE2 gene, which is a part of the planar cell polarity signalling pathway involved in kidney development. The top SNP associated with anaemia mapped to the intron of the putative PRKAG2 gene. The human ortholog of this gene has been associated with haematocrit and other blood-related traits, making it a prime candidate influencing parasite-triggered anaemia in brown trout. Our findings demonstrate the power of association mapping to pinpoint genomic regions and potential causative genes underlying climate change-driven parasitic disease resistance and severity. Furthermore, this work illustrates the first steps towards dissecting genotype-phenotype links in a wild fish population using closely related genome information.

Genomewide introgressive hybridization patterns in wild Atlantic salmon influenced by inadvertent gene flow from hatchery releases.

Many salmonid fish populations are threatened by genetic homogenization, primarily due to introgressive hybridization with hatchery-reared conspecifics. By applying genomewide analysis using two molecular marker types (1986 SNPs and 17 microsatellites), we assessed the genetic impacts of inadvertent gene flow via straying from hatchery releases on wild populations of Atlantic salmon in the Gulf of Finland, Baltic Sea, over 16 years (1996-2012). Both microsatellites and SNPs revealed congruent population genetic structuring, indicating that introgression changed the genetic make-up of wild populations by increasing genetic diversity and reducing genetic divergence. However, the degree of genetic introgression varied among studied populations, being higher in the eastern part and lower in the western part of Estonia, which most likely reflects the history of past stocking activities. Using kernel smoothing and permutation testing, we detected considerable heterogeneity in introgression patterns across the genome, with a large number of regions exhibiting nonrandom introgression widely dispersed across the genome. We also observed substantial variation in nonrandom introgression patterns within populations, as the majority of genomic regions showing elevated or reduced introgression were not consistently detected among temporal samples. This suggests that recombination, selection and stochastic processes may contribute to complex nonrandom introgression patterns. Our results suggest that (i) some genomic regions in Atlantic salmon are more vulnerable to introgressive hybridization, while others show greater resistance to unidirectional gene flow; and (ii) the hybridization of previously separated populations leads to complex and dynamic nonrandom introgression patterns that most likely have functional consequences for indigenous populations.

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin.

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of FST and to evaluate whether highly variable microsatellite loci resulted in deflated FST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (FCT = 0.010) and within (FST = 0.014-0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (PCT > FCT ). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (PST > FST ), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated FST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for PST -FST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.

Molecular identification of cryptic bumblebee species from degraded samples using PCR-RFLP approach.

The worldwide decline and local extinctions of bumblebees have raised a need for fast and accurate tools for species identification. Morphological characters are often not sufficient, and molecular methods have been increasingly used for reliable identification of bumblebee species. Molecular methods often require high-quality DNA which makes them less suitable for analysis of low-quality or older samples. We modified the PCR-RFLP protocol for an efficient and cost-effective identification of four bumblebee species in the subgenus Bombus s. str. (B. lucorum, B. terrestris, B. magnus and B. cryptarum). We used a short partial mitochondrial COI fragment (446 bp) and three diagnostic restriction enzymes (Hinf I, Hinc II and Hae III) to identify species from degraded DNA material. This approach allowed us to efficiently determine the correct species from all degraded DNA samples, while only a subset of samples 64.6% (31 of 48) resulted in successful amplification of a longer COI fragment (1064 bp) using the previously described method. This protocol can be applied for conservation and management of bumblebees within this subgenus and is especially useful for fast species identification from degraded samples.

Genetic mixed-stock analysis of lake-run brown trout Salmo trutta fishery catches in the Inari Basin, northern Finland: implications for conservation and management.

Genetic mixed-stock analysis (MSA) of wild lake-run brown trout Salmo trutta fishery catches (n = 665) from the Inari Basin (northern Finland) between 2006 and 2008 was carried out using a previously characterized baseline with 30 populations (n = 813) and 13 microsatellite loci. Altogether, 12 populations contributed significantly to mixed-stock fisheries, with the Ivalojoki system being the major contributor (70%) to the total catch. When catches were analysed regionally, geographically nearby populations were the main contributors to the local catches, indicating that a large proportion of S. trutta occupy lacustrine areas near the natal river mouth rather than dispersing throughout the lake. Similarly, far upstream populations contributed insignificantly to catches. These findings have important implications for the conservation and sustainable fishery management of the Inari system.

Permanent genetic resources added to Molecular Ecology Resources Database 1 February 2012 - 31 March 2012.

This article documents the addition of 171 microsatellite marker loci and 27 pairs of single nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Bombus pauloensis, Cephalorhynchus heavisidii, Cercospora sojina, Harpyhaliaetus coronatus, Hordeum vulgare, Lachnolaimus maximus, Oceanodroma monteiroi, Puccinia striiformis f. sp. tritici, Rhea americana, Salmo salar, Salmo trutta, Schistocephalus solidus, Sousa plumbea and Tursiops aduncus. These loci were cross-tested on the following species: Aquila heliaca, Bulweria bulwerii, Buteo buteo, Buteo swainsoni, Falco rusticolus, Haliaeetus albicilla, Halobaena caerulea, Hieraaetus fasciatus, Oceanodroma castro, Puccinia graminis f. sp. Tritici, Puccinia triticina, Rhea pennata and Schistocephalus pungitii. This article also documents the addition of 27 sequencing primer pairs for Puffinus baroli and Bulweria bulwerii and cross-testing of these loci in Oceanodroma castro, Pelagodroma marina, Pelecanoides georgicus, Pelecanoides urinatrix, Thalassarche chrysostoma and Thalassarche melanophrys.

High level of population genetic structuring in lake-run brown trout, Salmo trutta, of the Inari Basin, northern Finland.

Rivers draining into (Lake) Inarijärvi, northern Finland, sustain a number of lake-run brown trout, Salmo trutta, populations but, as with most lake-run S. trutta systems, the level of population genetic structuring among populations is unknown. To address this and to assist fish stock management in the region, the population genetic structure of S. trutta collected from 28 sampling sites in rivers flowing into Inarijärvi was studied using 13 microsatellite loci. Populations were clustered into three separate groups, largely corresponding to geographic regions, with between-region F(ST) values ranging from 0·11 to 0·16. The significant differentiation observed between most populations within each region also implies that individual populations should be recognized as separate management units and actions to improve, and subsequently maintain, conditions for natural spawning should be prioritized. The results of this study further indicate that the trout from each of these regions may have different biological characteristics, such as local-lake feeding behaviour among the western populations and strong isolation among the northern stocks. As a consequence, further research is warranted to better understand the level of ecological uniqueness of lake-run S. trutta populations.

Beyond MHC: signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci.

Using Atlantic salmon (Salmo salar) as a model system, we investigated whether 18 microsatellites tightly linked to immune-relevant genes have experienced different selection pressures than 76 loci with no obvious association with immune function. Immune-relevant loci were identified as outliers by two outlier tests significantly more often than nonimmune linked loci (22% vs. 1.6%). In addition, the allele frequencies of immune relevant markers were more often correlated with latitude and temperature. Combined, these results support the hypothesis that immune-relevant loci more frequently exhibit footprints of selection than other loci. They also indicate that the correlation between immune-relevant loci and latitude may be due to temperature-induced differences in pathogen-driven selection or some other environmental factor correlated with latitude.

High Gyrodactylus salaris infection rate in triploid Atlantic salmon Salmo salar.

We describe an unusually high infection rate of Gyrodactylus salaris Malmberg in juvenile Atlantic salmon Salmo salar L. of Baltic Sea origin, which are generally believed to be more resistant to G. salaris than East Atlantic salmon populations. Based on analyses of mitochondrial (complete cytochrome oxidase 1 [CO1] gene, 1548 bp) and nuclear (ADNAM1, 435 bp; internal transcribed spacer [ITS] rDNA region, 1232 bp) DNA fragments, the closest relatives of the characterized Estonian G. salaris strain were parasites found off the Swedish west coast and in Raasakka hatchery, Iijoki (Baltic Sea, Finland). Analyses of 14 microsatellite loci of the host S. salarrevealed that approximately 40% of studied fish were triploids. We subsequently identified triploid Atlantic salmon of Baltic origin as more susceptible to G. salaris infection than their diploid counterparts, possibly due to compromised complement-dependent immune pathways in triploid salmon. This is in accordance with earlier studies that have shown elevated susceptibility of triploids to various viral or bacterial pathogens, and represents one of the first reports of increased susceptibility of triploid salmonid fish to an ectoparasite. However, further experimental work is needed to determine whether triploid Atlantic salmon is generally more susceptible to G. salaris compared to their diploid counterparts, irrespective of the particular triploidization method and population of origin.

Use of differential expression data for identification of novel immune relevant expressed sequence tag-linked microsatellite markers in Atlantic salmon (Salmo salar L.).

Despite the large number of genes contributing to the immune response, wildlife immunogenetic studies have tended to focus mostly on the major histocompatibility complex-related genes. Here, we utilized previously published microarray and competitive RNA hybridization information to identify 3750 immune relevant Atlantic salmon (Salmo salar) expressed sequence tags. We then identified those expressed sequence tags containing microsatellites and subsequently designed 48 primer pairs and tested them for polymorphism in Atlantic salmon. Altogether, 16 polymorphic markers were characterized, with allele numbers ranging from two to 18, and these 16 loci were further tested in five other salmonid species.

Challenges for identifying functionally important genetic variation: the promise of combining complementary research strategies.

Strategies for the identification of functional genetic variation underlying phenotypic traits of ecological and evolutionary importance have received considerable attention in the literature recently. This paper aims to bring together and compare the relative strengths and limitations of various potentially useful research strategies for dissecting functionally important genetic variation in a wide range of organisms. We briefly explore the relative strengths and limitations of traditional and emerging approaches and evaluate their potential use in free-living populations. While it is likely that much of the progress in functional genetic analyses will rely on progress in traditional model species, it is clear that with prudent choices of methods and appropriate sampling designs, much headway can be also made in a diverse range of species. We suggest that combining research approaches targeting different functional and biological levels can potentially increase understanding the genetic basis of ecological and evolutionary processes both in model and non-model organisms.

Extensive immigration from compensatory hatchery releases into wild Atlantic salmon population in the Baltic sea: spatio-temporal analysis over 18 years.

Genetic homogenization has been recognized as a serious threat in an increasing number of species, including many salmonid fishes. We assessed the rate and impact of immigration from the main hatchery stocks of Atlantic salmon in the Gulf of Bothnia into one of the largest wild salmon populations in the Baltic Sea, the River Vindelälven, within a temporal framework of 18 years (from 1985-2003). We provide genetic evidence based on mtDNA and microsatellite markers, using mixed-stock analysis, that a large proportion (66%) of fin-damaged spawners (n=181) caught in the Ume/Vindelälven during 1997-2003 originated from the hatcheries in the Rivers Angermanälven, Luleälven and Ljusnan. The maximum-likelihood estimate of immigration rate from these hatcheries into the wild Vindelälven population was 0.068 (95% CI 0.021-0.128) over the studied time period (1985-2003) and reached up to a quarter (m=0.249, 95% CI 0.106-0.419) of the total population during 1993-2000. This resulted in significant (P<0.01) genetic homogenization trend between the wild Vindelälven population and hatchery stocks of the Angermanälven and Luleälven. Our results demonstrate extensive straying from geographically distant hatchery releases into wild salmon population and emphasize the genetic risks associated with current large-scale stocking practices in the Baltic Sea.

Matrilinear phylogeography of Atlantic salmon (Salmo salar L.) in Europe and postglacial colonization of the Baltic Sea area.

Sixty-four samples from 46 salmon populations totalling 2369 specimens were used for polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis of the mitochondrial ND1 region. The final analyses included 3095 specimens from 60 populations in Northern Europe. A subsample was analysed by RFLP of ND3/4/5/6. Representative RFLP haplotypes from different parts of the distribution area were sequenced and the phylogeny of European haplotypes and their relations to the North American lineage was described. The four common European haplotypes derive from the ancestral ND1-BBBA (rooting the European clade to the North American) by one-step substitutions: AAAA < AABA < BBBA > BBBB. The Swedish west-coast populations differ from the geographically close southern Baltic, indicating absence of inward and limited outward gene flow through the Danish straits during the last 8000 years. Within the Baltic Sea, only three ND1 haplotypes were detected and there was no variation for ND3/4/5/6. In the whole southern Baltic and in lakes Vänern, Ladoga and Onega the haplotype AABA dominated. Proposed postglacial colonization routes to the Baltic Sea are discussed in relation to the haplotype distribution pattern.

Identification of the origin of an Atlantic salmon (Salmo salar L.) population in a recently recolonized river in the Baltic Sea.

The founder event in a recently recolonized salmon population in the Baltic Sea (Gulf of Finland) was investigated. To identify the origin of the founders, four wild populations and two hatchery stocks were analysed using six microsatellite loci. The results of assignment tests and factorial correspondence analysis suggest that the initial recolonizers of the river Selja originated from the geographically nearest (7 km) wild population (river Kunda) but as the result of stocking activities, interbreeding between recolonizers and hatchery individuals has occurred in subsequent years. Although the hatchery releases are outnumbering the wild salmon recruitment in the Baltic Sea at present, our results suggest that the native populations may still have an important role in colonization processes of the former salmon rivers.